Density and Temperature of Bosons from Quantum Fluctuations 
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A method to determine the density and temperature of a system is proposed based on quantum 
fluctuations typical of Bosons in the limit where the reached temperature T is close to the critical 
temperature T c for a Bose condensate at a given density p. Quadrupole and particle multiplicity 
fluctuations relations are derived in terms of This method is valid for weakly interacting infinite 
and finite Boson systems. As an example, we apply it to heavy ion collisions using the Constrained 
Molecular Dynamics (CoMD) approach which includes the Fermi statistics. The model shows some 
clusterization into deuteron and a clusters which could suggest a Bose condensate. However, our ap- 
proach demonstrates that in the model there is no Bose condensate but it gives useful informations 
to be tested experimentally. We stress the differences with methods based on classical approxi- 
mations. The derived 'quantum' temperatures are systematically higher than the corresponding 
'classical' ones. The role of the Coulomb charge of fragments is discussed. 

PACS numbers: 25.70.Pq,42.50.Lc, 64.70.Tg 
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Fragmentation of heavy ions displays an anomalous 
production of a particles as compared to nucleons [l|- 
3] . This poses the question of what is the role of Bosons 
in nuclear matter and finite nuclei. We know that light 
nuclei display a a-cluster structure which could be ex- 
emplified by the so-called 'Hoyle' state in 12 C, i.e. the 
first excited state of such a nucleus which decays into 
3a's(3|. The fact that the ground state of nuclei could 
be made of a clusters could justify its copious produc- 
tion in heavy ion collisions near the Fermi energy. At 
the same time these facts arise the natural question if a 
clustering and production could be signatures of a Bose- 
Einstein condensate. In this work we discuss some prop- 
erties at finite temperatures assuming either a classical 
gas or a quantum Bose system, a Fermi system has been 
discussed inQ. We show that at the densities and tem- 
peratures of interest the classical approximation is not 
valid. This is at variance with many experimental and 
theoretical results in heavy ion collisions near the Fermi 
energy [H, which assume the classical approxima- 

tion to be valid. We base our approach on fluctuations 
estimated from an event by event determination of frag- 
ments arising after the energetic collision. A similar 
method has recently been applied to observe suppres- 
sion of fluctuations in a trapped Fermi gas [l2| and in 
Bose condensate s [13 . In ref.[f|, we proposed a method 
to go beyond ref. |l2l Il3| by including quadrupole fluctu- 
ations as well to have a direct measurement of densities 
and temperatures for subatomic systems for which it is 
difficult to obtain such informations in a direct way. We 
apply our proposal to the microscopic CoMD approach 
[14j which includes Fermionic statistics. Because of an- 
tisymmetrization, the model gives some clustering into 
a like structure in the ground state of some nuclei such 
as 40 Co. Also, in fragmentation reactions, a larger yield 
of a clusters is observed, however the experimental yield 
is largely underestimated[f4[ • These features should be 
kept in mind when discussing a possible Bose condensate 
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Figure 1: Normalized variance versus exitation energy per 
nucleon 5]. (Top panel) CoMD results for d (full circles) and 
a particles (full squares). For comparison the normalized fluc- 
tuations for fermions (bottom panel). (Open) Circles, squares 
and triangles refer to protons, neutrons and tritons, stars refer 
to 3 He. Notice the change of scales in the two panels. 



in the model. In fact we do not find any clear evidence 
of the condensate in the model but we can intuitively 
understand how to obtain it and especially how to pin 
down the condensate in experimental data. More refined 
models are possible but experimental data are needed in 
order to guide the modeling. We believe that such data 
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Figure 2: Reduced density versus reduced temperature for 
Bosons. Symbols as in figure f . 



could be obtained from heavy ion collisions using An de- 
tectors and performing a careful event by event analysis, 
the major serious problem we foresee is in the event se- 
lection for which the results discussed here in terms of 
the CoMD approach could be of guidance. 
A method for measuring the temperature was proposed 
in Q based on momentum fluctuations of detected par- 
ticles. A quadrupole Q xy =< p% — pt > is defined in a 
direction transverse to the beam axis (z-axis) to minimize 
non equilibium effects [|| and the average is performed, 
for a given particle type, over events. Such a quantity 
is zero in the center of mass of the equilibrated emitting 
source. Its variance is given by the simple formula: 
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where f(p) is the momentum distribution of particles. In 
@ a classical Maxwell-Boltzmann distribution of parti- 
cles at temperature T c i was assumed which gives: o~ 2 y = 
N4m 2 Tci, m is the mass of the fragment. N is the av- 
erage number of particles which could be conveniently 
normalized to one. In heavy ion collisions, the produced 
particles do not follow classical statistics thus the correct 
distribution function must be used in eq.(l). Protons(p), 
neutrons (n). tritium(t) etc. follow the Fermi statistics 
whileQ, deuterium(d), alpha(a) etc., even though they 
are constituted of nucleons, should follow the Bose statis- 
tics. In this work we will concentrate on Bosons only and 
in particular d and a which are abundantly produced in 
the collisions thus carrying important informations on 
the densities and temperatures reached. Using a Bosc- 



Einstein distribution f(p) and expanding near the critical 
temperature T c at a given density p, we get (l5| : 
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where the g n { z ) functions are well studied in the 
literature fl~5j| and z is the fugacity which depends on the 
critical temperature for Bose condensation and thus on 
the density of the system [ID]. Notice the similarity with 
the classical result which is modified by the ratio of the 
g n (z) functions only. Below the critical temperature such 

a ratio is a constant 97/2 ^l =????, thus the Bose temper- 

53/2(0) ' ^ 

ature is higher than the classical one. Above the critical 
temperature the ratio of the g functions approaches one 
for large T and small densities thus recovering the clas- 
sical result. These features should be contrasted with 
Fermion systems for which the temperature is always 
smaller than the classical limitQ, i.e. opposite to the 
Boson case. These results are very important and could 
be used to pin down a Bose condensate by comparing 
Fermions and Bosons produced in nuclear reactions on 
an event by event basis. The quadrupole fluctuations de- 
pend on temperature and density through T c , thus we 
need more informations in order to be able to determine 
both quantities for T > T c . 

Within the same framework we can calculate the fluc- 
tuations of the d, a multiplicity distributions. These are 
given by [l5j : 
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The difference with the Fermionic case [5] is striking: for 
Bosons, fluctuations are larger than the average and 
might diverge near the critical point, eq.(5), in the indi- 
cated approximations. Interactions and finite size effects 
will of course smoothen the divergence [15]. Two solutions 
are possible depending if we are above or below the criti- 
cal temperature for a Bose condensate. Below the critical 
point, eq.(2) can be used to calculate T and then eq.(4) 
gives the critical temperature and the corresponding den- 
sity. Above the critical point it is better to estimate the 
chemical potential which, in the same approximation, is 
given by: 
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Notice the similarity of this result to eq.(3) in where 
the chemical potential is substituted by the Fermi energy. 
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Figure 3: Excitation energy per particle divided by temper- 
ature vs temperature. Symbols as in figure 1, open symbols 
refer to calculations without Coulomb. 



From this equation we can estimate the g n functions en- 
tering eq. (3) and obtain the value of T. Using such a value 
in eq.(5), gives T c and the density p. In the numerical 
simulations discussed below we can always use the two 
solutions but one of these can be rejected using physical 
considerations. For instance assuming that we are below 
the critical point leads to densities as high as ten times 
the ground state density which is unphysical for heavy 
ion collisions around the Fermi energy. 

To illustrate the strength of our approach we simulated 
40 Ca + 40 Ca heavy ion collisions at fixed impact param- 
eter b = lfm and beam energies Ei^/A ranging from 4 
MeV/A up to 100 MeV/A. Collisions were followed up to 
a maximum time t = 1000/m/c in order to accumulate 
enough statistics. Because of the low yield of fragments, 
more than 200,000 events per beam energy were calcu- 
lated. In order to correct for collective effects as much as 
possible, we defined a 'thermal' energy as in ref.Q. 

In figure 1 we plot the reduced variances versus excita- 
tion energy per particle[5|- The CoMD results are given 
by the full symbols, top panel. As we see in the fig- 
ure, a normalized fluctuations are generally larger than 
d-fluctuations. As we will show below, this implies that 
those particles might explore different regions of densities 
and temperatures. In both cases, fluctuations are large 
and, in some cases, above Poissonian for cv's. In order 
to understand if a Bose condensate occurs in the model 
(and in the future in experiments) it is instructive to 
compare the Bosons normalized fluctuations to those of 
Fermions discussed in ref.[5j. In figure 1-bottom panel, 
normalized Fermions fluctuations are given. As we see 



the normalized fluctuations of p and n are much smaller 
than 1 at variance with the Bosons case, which would 
suggest a condensate. However, heavier Fermion clus- 
ters such as 3 He and tritons, display fluctuations larger 
than d and smaller than a. In particular 3 He fluctua- 
tions are smaller than t which demonstrates that it is 
not 'directly' a Coulomb effect. These facts are impor- 
tant to understand what is happening in the model and 
eventually search for an experimental confirmation. We 
offer here an intuitive explanation of the relative role of 
normalized fluctuations for different particles type. The 
CoMD model is essentially classical with a constraint in 
the equation of motion which keeps the occupation prob- 
ability f(r,p,t) smaller than 1 as dictated by the Pauli 
principle for Fermions [H I2L Il4j). A further implementa- 
tion of the Pauli principle is in the collision term which 
avoids that colliding nucleons occupy phase space regions 
which are occupied by other nucleons. Thus the Pauli 
principle reduces the available phase space and in turn 
the normalized fluctuations. For this reason p and n fluc- 
tuations are smaller than Poissonian. When composite 
fragments are formed, d,t, etc., the effect of Pauli block- 
ing is reduced (also because those particles form at low 
densities, see below), thus fluctuations become compa- 
rable to their average value. The effect that reduces the 
available phase space is now the binding energy. Not all 
nucleons can form a bound state, especially if their rela- 
tive kinetic energies are larger than the potential one. For 
this reason d fluctuations are smaller than 3 He which are 
smaller than t and smaller than a fluctuations. In par- 
ticular, since the d is over bound in the model (about 7 
MeV) we expect the data to exhibit smaller fluctuations. 
On the other hand the binding energy of a is smaller than 
the data (about 20MeV) thus suggesting the possibility of 
larger fluctuations than displayed in figure 1. However, 
the important ingredient which is missing in the model is 
the possibility of boson-boson collisions ( a- a, d-d, etc.) 
and correlations. If such collisions occur their probabili- 
ties will be enhanced by the Bose factors (1 + f(r,p,t)) 
in contrast to the Pauli blocking factors 0, This will 
produce fluctuations larger than Poissonian, which is a 
signature of a Bose condensate. Of course the last word 
is in the experimental data, similar to trapped Fermi and 
Boson gases [HI EH- Notice in figure 1 the occurrence of 
a minimum at a similar excitation energy for p,n and d 
but not for heavier clusters. This will have an effect on 
the EOS as we will show below. 

It is interesting to discuss the densities 'seen' by the 
different Bosons during the reaction. A plot of density 
(divided by the ground state density) versus tempera- 
ture (divided by the critical temperature for a conden- 
sate) is given in figure 2. Notice the peculiar behavior 
of d and a clusters. While the latter are formed at a 
constant reduced temperature but at different densities 
for each beam energy, the deuterons are formed always 
at a very small (constant) density but at different tem- 
peratures. As we noticed above two effects are at play. 
The first is that there is no Pauli blocking for nucleons 
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Figure 4: Energy density versus temperature, top panel with 
Coulomb, bottom panel without Coulomb. Symbols as in 
figure 3. 



inside the clusters, the second is the different binding en- 
ergies. Since d are over bound in the model we expect to 
see that even smaller densities will be seen in the data, 
the opposite we expect for a. These features remind of 
Mott transitions and in particular suggest that different 
particle types might be sensitive to different regions of 
the NEOS as already noticed for Fermion -V. This is an 
interesting results since it gives the opportunity to de- 
termine the EOS even though we are dealing with small, 
dynamical systems, where we need to correct for some 
effects such as Coulomb. 

In figure 3 we plot the thermal energy divided by T versus 



temperature. From eq.(2-3) we expect to get the classical 
result at very high T and this is the reason for the y-axis 
choice in the figure. As we see the results indeed approach 
the value of 3/2 as expected, even though we notice a de- 
crease at the largest temperature which is an indication 
of transparency in the model. For some calculations we 
have also turned off the Coulomb field (open symbols). 
It is evident that the general features are the same apart 
from a shift in T, similar to what was observed for p and 
nQ. It is interesting to notice that a particles indicate 
a first order phase transition, while d might suggest a 
crossover or a second order phase transition. Recall that 
the densities explored by the two fragments are quite 
different, thus it is not surprising if they behave differ- 
ently. These might be better grasped in figure 4 where 
the energy density is plotted versus temperature. This 
figure should be compared to the Fermion caseQ where a 
rapid increase is noticed at about 3-4 MeV temperature. 
Such an increase is seen in the d-case but there is a deep 
minimum near such temperatures, reflecting the fluctua- 
tions displayed in figure 1. As we see in figure 4 (bottom 
panel) , such a behavior is not due to Coulomb which sim- 
ply shifts the temperatures to higher values. Of course 
a mixture of bosons and fermions produces a quite rich 
EOS as in the case of 3 He, 4 He gas mixtures[15], thus 
our proposed approach might open a quite interesting 
field of research. 

In conclusion, in this work we have addressed a gen- 
eral approach for deriving densities and temperatures of 
bosons. In the framework of the Constrained Molecular 
Dynamics model, we have discussed collisions of heavy 
ions below lOOMeV/A and obtained densities and tem- 
peratures at each bombarding energy. Knowing the ther- 
mal energy of the system, we can derive the energy den- 
sity and temperature reached during the collision. We 
have seen in this work that different particles like d and 
a explore different density and temperature regions, sim- 
ilar to the Fermion case[5|. Open problems such as a 
Bose condensate in nuclei, Mott transitions, pairing etc. 
in low density matter might be addressed through a de- 
tailed study of the EOS. 
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